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INTRODUCTION 

The Mariner 4 S-band occultation experiment (Kliore et al., 1965) and ground-based spectro- 
scopic observations (Owen, 1966; Spinrad et al., 1966; Belton and Hunten, 1966) established that the 
pressure at the surface of Mars is typically less than 10 mb and that C0 2 is the major constituent of 
the atmosphere. From its thinness and composition, one can conclude that the Martian atmosphere 
approaches radiative equilibrium at low altitudes, with diurnal temperature adjustments near the sur- 
face (Gierasch and Goody, 1968). The results of the preliminary analysis of the Mariner 6 and 7 
occultation data for the lower atmosphere (Kliore et al., 1969) were largely inconsistent with this 
conclusion, however, with several features of the derived temperature profiles differing surprisingly 
from theoretical predictions. These results seemed to indicate that important and hitherto unsus- 
pected physical processes might be at work in the lower atmosphere of Mars. For example, at the 
Mariner 6 entry point, an unexpectedly warm, nearly isothermal layer underlying a region of precipi- 
tously falling temperature was found in the middle atmosphere. Neither of these features had ever 
been produced in any thermal model, and to account for their presence, considerable modification of 
the widely accepted radiative-convective theory of the lower Martian atmosphere would be required. 

In an attempt to resolve this discrepancy between theory and observation, which has an impor- 
tant bearing on our understanding of the meteorology of Mars, an independent analysis of the Mariner 
6 and 7 occultation data has been performed at Goddard Space Flight Center’s Institute for Space 
Studies (GISS). A reanalysis of the data was also conducted at the Jet Propulsion Laboratory (JPL) 
by the experimental team, with results described by Kliore et al. (1970). 



The GISS investigation has involved the construction of two separate computer programs for— 

(1) Combining Doppler residual data with trajectory data to obtain the raw residuals as a func- 
tion of the distance of the Mariner spacecraft from the limb of Mars. 

(2) Filtering from the raw residual data those components not due to the Martian atmosphere. 

(3) Integrating the filtered residuals over time to derive the phase change as a function of the 
spacecraft’s distance from the planetary limb. 

(4) Inverting the phase-change data, with a suitable correction for refractive bending, to obtain 
profiles of refractivity in the Martian atmosphere. 

(5) Deriving temperature, pressure, and electron-density profiles on Mars from the refractivity 
distributions. 

(6) Computing the systematic and formal errors introduced at each of the above stages of the 
analysis. 

Earlier (Rasool et al., 1970), the authors described their partial and preliminary results concern- 
ing the temperature distributions in the lower Martian atmosphere at the four occultation points. 
Those results differed in several important respects from the preliminary findings of Kliore et al. 
(1969). For example, the warm, nearly isothermal region in the middle atmosphere found in the 
Mariner 6 entry temperature profile of Kliore et al. (1969) was not present in the corresponding pro- 
file of Rasool et al. (1970). However, agreement with Kliore et al. (1969) was obtained with regard 
to the suggestion of extremely low temperatures at higher altitudes on Mars. Also, the results of 
Rasool et al. (1970) were in general agreement with the set of corrected thermal profiles presented by 
Kliore et al. (1970). Discrepancies between the first results of Kliore et al. (1969) and the later find- 
ings of the JPL and GISS groups were attributed to unsatisfactory correction of the Doppler data in 
the preliminary analysis. 

The authors have now obtained final profiles of temperature, pressure, and electron density on 
Mars from the occultation data for the Mariner 6 and 7 entry and exit cases, and the results are pre- 
sented for both the lower atmosphere and the ionosphere of Mars. Also included are the results of a 
detailed analysis of the systematic and formal errors introduced at each stage of the data-reduction 
process. 

GENERAL PRINCIPLES AND REQUIRED ASSUMPTIONS 

The possibility of determining certain properties of a planetary atmosphere by monitoring the 
telemetry link between the Earth and a space probe on an occulting trajectory was first discussed by 
Fjeldbo and Eshleman (1965). As the probe approaches the limb of the planet, the phase path, fre- 
quency, and amplitude of the signal are altered by its passage through the planetary atmosphere. The 
local index of refraction differs from unity in the neutral atmosphere and ionosphere, which causes 
the phase and group velocities to depart from their free-space values. As the mean refractive index 
along the ray path varies with the distance of the path from the limb of the planet, differential bend- 
ing of the rays occurs, producing focusing and defocusing of the signal. 
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The changes in the signal characteristics are a function of the distribution of the index of refrac- 
tion in the atmosphere of the planet. In principle, if these changes are recorded as the spacecraft 
passes behind the planet, the distribution of refractive index can be inferred. 

The present investigation is concerned with the determination of atmospheric parameters on Mars 
from the Mariner 6 and 7 phase-change data. The phase changes produced by the Martian atmosphere 
must be obtained by subtracting all predictable phase shifts due to light transit time, spacecraft motion, 
Earth rotation, and so forth from the total phase change during the period of the experiment. Since 
the effect of the atmosphere is small in comparison with these other effects, high accuracy is required 
at this stage in the data-reduction process; of prime importance is a precise knowledge of the space- 
craft trajectory and range rate. If the necessary accuracy is achieved, model fitting or, alternately, 
integral inversion of the residual phase changes can be used to estimate the properties of the atmos- 
phere. 

The residual phase change is in general a result of the combined effects of the neutral atmosphere 
and ionosphere of the planet. Since the contribution of the neutral atmosphere to the phase change is 
proportional to the total atmospheric density along the ray path, this component will be largest when 
the ray path is close to the limb of the planet. On the other hand, the effects of the ionosphere on 
the phase change are proportional to the integrated ion content along the ray path, and these effects 
will be greatest when the ray path lies at a distance from the limb approximately equal to the height 
of the ionosphere above the surface. If the neutral and ionized regions which produce a detectable 
change in the signal characteristics are separated in altitude, it should be possible to obtain informa- 
tion on both regions from a single-frequency experiment. 

Obviously, the measured variation of residual phase change with distance from the planetary limb 
does not uniquely determine the spatial distribution of refractive index in the atmosphere; conse- 
quently, some a priori assumptions must be made in interpreting the atmospheric phase-change data. 
The assumption of a spherically symmetric, stratified atmosphere simplifies model fitting and permits 
the use of integral inversion techniques in deriving a refractivity distribution from the observed phase 
changes. Since a highly localized region of the lower atmosphere is probed by the telemetry signal, 
the assumption of spherical symmetry is equivalent, in this case, to postulating that vertical gradients 
of refractivity in the lower atmosphere are much larger than horizontal gradients— an assumption 
often used in atmospheric physics. The assumption of spherical symmetry in the Martian ionosphere 
is difficult to justify, however, and the effects of possible departures from this condition in the region 
probed on the interpretation of the phase data should be considered. For example, in the ionosphere, 
the electron densities are likely to be symmetric about the subsolar point; thus, the degree to which 
the nonspherical nature of the ionosphere affects the interpretation of the phase data should depend 
on the local zenith angle at the point of occultation. 

From the above discussion, it should be clear that only if the following three conditions are met 
can the single-frequency occultation experiment succeed in providing accurate information on the dis- 
tributions of temperature, neutral density, and electron concentration in the Martian atmosphere: 

(1) The changes in the telemetry signal that result from the Martian atmosphere can be isolated from 
the changes that result from all other effects; (2) a physical separation in altitude exists between the 
neutral and ionized regions which produce these changes; and (3) the distribution of refractivity in 
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the atmosphere must not contradict a basic assumption about the atmospheric structure, that horizon- 
tal gradients of refractivity are much smaller than vertical gradients. 

CORRECTION OF DOPPLER RESIDUAL INPUT DATA 

The present analysis is based upon the Mariner 6 and 7 closed-loop Doppler residual and trajec- 
tory data supplied by the JPL experimental team. During the near encounter of Mariner 6 with Mars 
at 0519 GMT, July 31, 1969, Deep Space Network (DSN) Stations 12 and 14 at Goldstone, California, 
monitored the spacecraft telemetry. At the time of the closest approach of Mariner 7 (0501 GMT on 
August 5, 1969), Station 41 at Woomera, Australia, was also in position to communicate with that 
spacecraft, and all three stations monitored the telemetry in this case. Ten distinct sets of Doppler 
residual data are available, including two each (Stations 12 and 14) for Mariner 6 occultation entry and 
exit, and three each (Stations 12, 14, and 41) for Mariner 7 occultation entry and exit. 

Four additional bodies of data with slightly reduced noise levels can be obtained by averaging the 
data from the different stations in operation during the Mariner 6 and 7 encounters. This averaging 
procedure reduces the level of the noise that originates independently at each receiving station during 
the experiment, such as the effects of changes, due to turbulence, in the refractivity along the ray path 
through the Earth’s atmosphere. Averaging over stations does not remove or reduce the component of 
the noise produced by the spacecraft systems, which must be identical at all stations. 

Entry data were taken in a two-way mode, with a frequency generated by a rubidium standard 
transmitted to the spacecraft, which coherently retransmitted it to Earth. Exit data were acquired in 
a one-way mode, with a crystal oscillator onboard the spacecraft providing the reference signal for 
transmission to Earth. As a consequence of the two different modes of operation, immersion data 
reflect a two-way passage through the atmosphere of Mars, whereas emersion data reflect a one-way 
passage. Also, according to Fjeldbo et al. (1970), the stability of the rubidium standard used in acquir- 
ing the entry data (a few parts in 1 0 1 2 ) is greater by nearly two orders of magnitude than the stability 
of the onboard crystal oscillators used during exits (about one part in 10 10 ). 

Both the Doppler residual data points and the trajectory points have been supplied at 1-s time 
intervals. However, the trajectory data are shifted in time from the residual data by a fraction of a 
second; thus, the position of the spacecraft at the times at which the residual points are given must be 
obtained by interpolation. This interpolation was performed in the present analysis through the use 
of a third-order-polynomial spline fit to the trajectory data. 

Also, in each case, it was necessary to extrapolate the Doppler residual data to obtain the value 
of the residual at the precise time of occultation entry or exit. Since both spacecraft approached 
occultation with velocities close to 3 km/s perpendicular to the Earth-Mars line, a careful extrapolation 
of the residuals is of the utmost importance for an accurate determination of the atmospheric condi- 
tions at the surface of Mars. Because the refractivity decreases exponentially with altitude above the 
surface of the planet, it can be anticipated that a rapid variation of the residuals will occur immediately 
prior to and after occultation. To obtain a value for the Doppler residual at the moment of occulta- 
tion, the 10 residual points preceding or following occultation were least-squares fitted with a third- 
order polynomial. This polynomial was extrapolated (over an interval of less than 1 s) to the occulta- 
tion time to obtain the residual value at the instant of occultation. 
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The input Doppler residual data are characterized by large quasi-linear drifts during the periods 
before and after occultation of the spacecraft by Mars. Almost identical patterns are found in the two 
or three different sets of data corresponding to a single entry or exit point. As an example, Figures 1 
through 4 show the raw Doppler residual data from Station 1 2 for the four occultation points on Mars. 
Except for minute differences due to random noise (which can presumably be eliminated by the 
averaging described above), the data from the other stations are indistinguishable to the eye from the 
Station 1 2 data shown in these figures. 

The trends evident in the input data of Figures 1 through 4 are probably due to the combined 
effects of oscillator drifts and trajectory errors.* Because they extend far into preencounter and post- 
encounter periods, it is clear that they are in no way related to the atmospheric characteristics on Mars. 
However, it is not possible to point to the exact cause of these drifts in the residual data; thus, their 
functional form cannot be specified with certainty. 

To determine the properties of the Martian atmosphere from the residual data, the unwanted 
effects evident in Figures 1 through 4 were removed by subtraction of least-squares-fitted bias curves. 

In our preliminary analysis, a fitted curve of the form A + Bt + C sin t + D cos t was subtracted from 
the data (Rasool et al., 1970). Two additional avenues of approach have been used in the present study. 
In one test, Fourier series of varying complexity were fitted to the data in the preencounter and post- 
encounter periods; in another, polynomials of different orders were fitted to the preencounter and 
postencounter data. 

For both Mariners 6 and 7, preencounter and postencounter periods were defined as the periods 
ending 130 s before occultation entry and beginning 40 s after occultation exit, respectively. With 
these definitions, the preencounter and postencounter periods end and begin, respectively, when the 
ray path is about 360 km from the limb of Mars along the perpendicular to the Earth-Mars line. 

Seven-term Fourier series and third-order polynomials can be fitted to all of the preencounter or 
postencounter Doppler-data sets with standard deviations of less than 0.1 Hz. This same degree of fit 
was found in the analysis reported earlier (Rasool et al., 1970). Although series with as many as 19 
terms and polynomials of as many as nine orders were considered in these experiments, no significant 
improvement of the fits to the preencounter and postencounter data was realized by considering series 
with more than seven terms or polynomials of order higher than third. The use of 19-term Fourier 
series and ninth-order polynomials produces a decrease of less than 0.02 Hz in the above standard 
deviations. 

Interestingly, the relative stability of the entry and exit frequency references described by 
Fjeldbo et al. (1970) is not reflected in the noise level of the data about the fitted curves. 

Unfortunately, difficulties arise when an attempt is made to extrapolate the higher-order fitted 
curves from the preencounter and postencounter periods into the regions where the atmospheric 
effects are present. In the case of the more complex functions considered, the addition of a single 
point to the group of data fitted can produce drastic changes in the extrapolated portion of the curve. 
This sensitivity is due to improper functional representation of the mechanisms that produce the 
drifts in the raw residual data. It is not difficult to imagine that the improper selection of one of these 


*Kliore, private communication. 
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Figure 2-Raw Doppler residual data (Station 12) for Mariner 6 
exit case. The abscissa is an arbitrary time scale. 
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Figure 3— Raw Doppler residual data (Station 12) for Mariner 7 
entry case. The abscissa is an arbitrary time scale. 
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Figure 4— Raw Doppler residual data (Station 12) for Mariner 7 
exit case. The abscissa is an arbitrary time scale. 



curves would result in a complete misrepresentation of the component of the raw Doppler data that is 
due to the Martian atmosphere. 

Changes in the extrapolated portions of the fitted functions as the sample of points fitted is varied 
are significantly smaller when simpler functions are considered. Nevertheless, on the basis of the some- 
what inconsistent behavior exhibited by even the simpler of the fitted functions, the problem of objec- 
tively removing the drifts from the raw residual data without having precise information on the mech- 
anisms that gave rise to these drifts would seem hopeless. 

To ensure the extraction of significant information on the Martian atmosphere from these bodies 
of data, it is necessary to apply an additional constraint on the acceptability of the fitted curves, 
beyond that of a small standard deviation in the preencounter or postencounter period. Such a 
criterion is suggested by one of the basic assumptions required for the success of the occultation 
experiment— that somewhere in the middle atmosphere of Mars, a region of essentially zero (undetecta- 
ble) refractivity exists. Unless this assumption is valid, the effects of the ionosphere and neutral atmos- 
phere on Mars cannot be unambiguously separated. 

With this assumption, a curve that properly represents the nonatmospheric component in the raw 
residual data must be ( 1 ) one that can be fitted with a small standard deviation to the data in the pre- 
encounter or postencounter region and (2) one whose removal from the raw residual data will even- 
tually produce a profile of refractivity on Mars, characterized by an undetectably low refractivity value 
over an extensive region of the middle atmosphere. This desirable feature of a refractivity profile will 
henceforth be referred to as a zero-refractivity ledge. (In this experiment, an undetectably low value 
of refractivity corresponds to an absolute value less than 0.05 «-unit). To derive curves with both of 
the above characteristics, the following procedure was followed. 

Polynomials of various orders were fitted to a particular group of preencounter or postencounter 
raw Doppler residuals. Each curve fitted was subsequently extrapolated into the region immediately 
prior to or after occultation where atmospheric effects are manifest and subtracted from the entire 
body of raw data to obtain a set of corrected residuals. These corrected residuals were then processed 
(as described below) to derive a distribution of refractivity as a function of the distance from the 
center of the planet to the point of closest approach of a ray. The refractivity profile was then exam- 
ined for consistency with the assumption of a zero-refractivity ledge. The ledge value for each case was 
determined by scanning for the minimum of the absolute value of the refractivity in the middle 
atmosphere. 

If consistency with the zero-ledge assumption for the Martian mesosphere was not found (refrac- 
tivity ledge of absolute value greater than 0.05 «-unit), the fitted curve was rejected as an unaccepta- 
ble representation of the nonatmospheric component of the raw residuals in the extrapolated region. 

If consistency was obtained (refractivity ledge of absolute value less than 0.05 w-unit), the refractivity 
distribution was used to derive distributions of temperature, pressure, and electron density for a chosen 
atmospheric composition. The entire procedure was then repeated, with another point being added to 
the group of preencounter or postencounter raw residuals to be fitted. 

A similar sequence of operations was used to select acceptable Fourier representations of the non- 
atmospheric component of the Doppler residual input data. 
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It should be noted that for all 1 0 individual sets of raw data and for all four sets of averaged data, 
several hundred different bias curves consistent with the zero-ledge assumption could be closely fitted 
to preencounter or postencounter data groups. All of the alternatives in each case produced closely 
similar results (described below) when the analysis was carried through to obtain profiles of tempera- 
ture, pressure, and electron density. 

Typical sets of residual data for the four points of occultation on Mars, corrected via the method 
outlined above, are shown in Figures 5 through 8. These correspond to averages over the stations that 
monitored the telemetry at each point, and they closely resemble the mean of the several hundred 
sets of corrected residual data that were found to produce refractivity profiles satisfying the zero-ledge 
criterion. 

DATA REDUCTION 

The basic geometry of the occultation experiment is illustrated in Figure 9. Fjeldbo and Eshleman 
(1968) have shown that this simple geometry approximates to a high degree of accuracy the actual 
occultation configuration for a thin atmosphere. At the time of entry into and exit from occultation, 
the Mariner 6 spacecraft was approximately 9700 km and 17 700 km behind Mars, respectively. In 
the case of Mariner 7, the distances behind the planet were 8500 km at immersion and 19 800 km at 
emersion. The Earth-Mars distances were about 95.7 million km and 100.2 million km, respectively, at 
the times of the Mariner 6 and 7 encounters. 

The analysis begins with a set of TV Doppler residuals 8., corrected by the procedure detailed 
above, and a set of TV spacecraft coordinates X Sj , Z Sj , corresponding to TV points of time t.. The index 
i increases and decreases with time in the cases of entry and exit, respectively; thus, point 1 corre- 
sponds to the farthest distance from the planet, the point at which the analysis is begun, and point TV 
corresponds to the surface of the planet. With reference to Figure 9, the altitude h t is given by 

<•> 

where X. is the distance from the center of Mars of the point of closest approach of a ray correspond- 
ing to the X s . f Z s . position of the spacecraft, and X p is the planetary radius at the occulation point. 

The phase change <j>. corresponding to the altitude of closest approach h i at time t. is obtained by 
time integration of the Doppler residuals according to the trapezoidal rule: 

/5i+5 2 \ (8 2 + 5 3 \ (8._ i+sA 

4>i = y—j— ) + \~r~ l j At 2 >3 + ■ " + \ — t~) (for l > 1} * (2) 

where . is the time interval between the Doppler residual measurements at times t._^ and t f . In 
the entry cases, integration was begun approximately 130 s before immersion, at the end of the pre- 
encounter period defined above, and carried forward in time. In the exit cases, the time order of the 
integration was reversed, with the integration starting at the beginning of the postencounter period, 

40 s after emersion, and extending backward into the period where atmospheric effects are present. 

In the bent-ray approximation of Figure 9, the distance of the closest approach of a ray from 
the center of the planet is related to the spacecraft coordinates and phase-path increase by 
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Figure 5— Raw Doppler residuals (dots), polynomial 
least-squares fit (smooth curve), and corrected resid- 
uals (noisy curve) corresponding to Mariner 6 entry. 
The residuals are averages of Station 12 and 14 data. 
The abscissa is an arbitrary time scale. 


Figure 6-Raw Doppler residuals (dots), polynomial 
least-squares fit (smooth curve), and corrected resid- 
uals (noisy curve) corresponding to Mariner 6 exit. 
The residuals are averages of Station 12 and 14 data. 
The abscissa is an arbitrary time scale. 
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Figure 8— Raw Doppler residuals (dots), polynomial 
least-squares fit (smooth curve), and corrected resid- 
uals (noisy curve) corresponding to Mariner 7 exit. 
The residuals are averages of Station 12, 14, and 
41 data. The abscissa is an arbitrary time scale. 
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Xj — X s . Z s . tan ct. , (3) 

where a f is the angle by which the path of the signal is bent from a straight line by the atmosphere: 


A5 f 


(4) 


where V Si is the component of the spacecraft velocity along the perpendicular to the Earth-Mars line 
and A is the wavelength of the telemetry signal (about 13.6 cm for both spacecraft). 

The straight-line approximation to the phase change $ is, with reference to Figure 9, 


<f>. = 0. — 




A cos a. 



(5) 


For small a., expanding cos a. and neglecting higher-order terms causes Equation 5 to reduce to 


<h. = 0. — 

i “i 



( 6 ) 


The straight-line portion of the total phase change is due to a region of nonzero refractivity, the Mar- 
tian atmosphere, encountered along the path of the telemetry signal. The contribution to the total 
phase change represented by the second term on the right-hand side of Equation 6, is associated with 
the bending of the propagation path from a straight line. 

The relationship of the straight-line phase change to the refractivity distribution in the atmosphere 
is given by 


10 -6 r + ~ 

*,= — / Rir)dZ, (7) 


where R(r) represents the refractivity of an element of unit volume located on the zth ray path at a dis- 
tance r from the center of the planet: 


r 2 = Xf +Z 2 . (8) 

Fjeldbo and Eshleman (1968) have outlined a numerical technique for obtaining the refractivity 
distribution in the case of a spherically symmetric Martian atmosphere, where R(r) is a constant. The 
geometry applicable in this case is illustrated in Figure 10, adapted from their Figure 18. Under this 
assumption, N rays pass through an atmosphere consisting of N spherical layers, each of constant 
thickness and constant refractivity. The integral of Equation 7 may then be replaced by a linear equa- 
tion of the form 
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Figure 9— Ray path geometry (adapted from Fjeldbo 
and Eshleman, 1965). 


Figure 10— Geometry applicable to a spherically 
symmetric Martian atmosphere (adapted from 
Fjeldbo and Eshleman, 1968). 


where R . is the refractivity of the / th layer and A Z. . designates one-half the distance traversed by the 
ith ray in the /th layer. Referring to Figure 10, one sees that the relationship of A Z. . to the radii of 
the atmospheric layers is 



when i is greater than /, and 


( 10 ) 
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The set of N linear equations of the form of Equation 9 may be solved for i? to yield the dis- 
tribution of refractivity in the atmosphere of Mars: 

R x =(X- 10 6 4 > 1 )/AZj ^ 


and 


7?, = ( X • 10 6 <i>. - 2y fl.AZ. -4- (for i > 2) . (12) 

\ h / “ 

To extract the distribution of neutral-particle number density in the lower atmosphere from this 
refractivity distribution, the chemical composition of the atmosphere must be specified. Ground-based 
infrared spectroscopic observations (Owen, 1966; Spinrad et al., 1966; Belton and Hunten, 1966) indi- 
cate that C0 2 , with a partial pressure of 5 ± 2 mb, is the major constituent of the Martian atmosphere. 
Such measurements also show that the actual and possible abundances for other gases which possess 
infrared spectra are too low for these gases to contribute significantly to the atmospheric refractivity. 
However, gases inactive in the infrared, such as N 2 , Ar, or Ne, could be present on Mars in significant 
amounts but would not be observable from the Earth’s surface. For example, according to the results 
of the ultraviolet spectroscopy experiments on Mariners 6 and 7, as much as 5 percent of the Martian 
atmosphere could be N 2 (Barth et al., 1969; Dalgamo and McElroy, 1970). 

It is assumed here that no gases other than C0 2 , N 2 , and Ne contribute to the observed refractiv- 
ity on Mars. The molecular refractivities n for the nonpolar molecules C0 2 and N 2 are 1.84 X 1 0 -1 7 
and 1.10 X 10 -17 , respectively, according to the data of Essen and Froome (1951). Tyler and Howard 
(1969) give a value similar to the above for the refractivity of C0 2 under Martian conditions. A value 
of 2.44 X 10 -18 for the molecular refractivity of Ne was derived from the refractive indexes tabulated 
in Astrophysical Quantities (1963). After specifying the fractional concentration /of each of these 
three atmospheric constituents at each level, the distribution of total number density N in the Martian 
atmosphere can be obtained from relationships such as 

Nj = R i l{nco 2 fco 2 . 2 /n 2 j + ”Ne/Ne,-) • (13) 


From the vertical distribution of total number density, the thermal structure of the Martian 
atmosphere can be obtained by integration of the hydrostatic law. The atmospheric temperature distri- 
bution in terms of the distribution of N is given by 


N(h 0 ) 1 f h o 

T(h) = T(h 0 ) — — + / gctmom) d$, (14) 

u N(h) kN(h) J 
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where h Q is the level selected as the upper boundary of the lower atmosphere, k is Boltzmann’s 
constant, g is the acceleration due to gravity, and in is the mean molecular mass. The level h Q is 
chosen to correspond to the position of the zero-refractivity ledge (discussed in the previous section). 

In the case of a homogenous atmosphere, the denominator of the right-hand side of Equation 1 3 
is a constant, and in is independent of altitude. Thus, Equation 1 4 simplifies to 

R(h 0 ) in r h o 

Tm - T ^-rn + im] i *«>*«><«■ <'« 

By using the trapezoidal rule to integrate the refractivities in Equation 1 5, we obtain (for i > / Q ) 

R >0 ™ & 

r r r -oI + ^^ /+ i^ /+ i)A/z ; , /+1 , (16) 

* l /'= f o 

where i and / Q are the indexes corresponding to altitudes h and h Q , respectively, and Ah. . +l is the 
height increment between levels / and / + 1 . 

The first term on the right-hand side of Equation 16 depends upon 7} 0 , the arbitrarily chosen 
temperature at the upper boundary of the lower atmosphere. This choice of boundary temperature, 
therefore, will to some degree influence the derived profile at lower altitudes. Because of its inverse 
dependence on the refractivity, however, the magnitude of this term decreases quite rapidly as deeper 
atmospheric levels are considered. The relationship between the refractivity distribution and the dis- 
tribution of free electrons in the upper atmosphere (the region above the zero-refractivity ledge) is 
given by 


Rv 2 

N ei = (cm' 3 ), (17) 

40.3 X10 1 2 

where v is the frequency of the telemetry signal (approximately 2.2 X 10 9 Hz). 

ERROR ANALYSIS 

Both systematic and random errors are involved in the derivation of the physical characteristics of 
the Martian atmosphere from the occultation data. Errors of the former type arise from oscillator drift 
and imperfect knowledge of the spacecraft trajectory. Errors of the latter kind result from the noise 
level of the input data. 

Our analysis assumes that the systematic errors are manifest in the large departures of the Doppler 
residuals from zero during the preencounter and postencounter periods, when the atmosphere on Mars 
cannot be modifying the spacecraft telemetry signal. These systematic errors are assumed to be 
removed, or largely reduced, by (1) performing a polynomial regression analysis on the data over 
selected time periods in the preencounter and postencounter phases, (2) extrapolating the resulting 
least-squares-fitted bias curve into the encounter period, and (3) subtracting this curve from the com- 
plete set of Doppler residuals to obtain a complete set of corrected residuals. 


15 



Since the error involved in extrapolating the bias curve through the encounter period cannot be 
formally estimated, we have placed the additional constraint of a near-zero-refractivity ledge, dis- 
cussed previously, on the data-reduction process. However, even this constraint does not guarantee 
that the extrapolation error is negligible; the authors have therefore analyzed a large number of cases 
resulting from different polynomial fits to the raw residuals to assess the validity of this correction 
procedure and to establish the range of the systematic errors that persist if the assumption of a low- 
refractivity region is correct. 

The polynomial regression analysis described above yields the variance of the Doppler residuals to 
which the bias curve is fitted. This variance is a measure of the noise level, or random error, inherent 
in the input data. As mathematical operations are performed, this initial random error will propagate 
through the analysis to all derived quantities. In contrast to the situation with regard to the systematic 
error described above, a formal procedure is available for calculating the resultant errors in quantities 
derived from the Doppler residuals. 

Our analysis of the formal errors accumulated during the derivation of the atmospheric charac- 
teristics is based on the following assumptions: 

(1) Systematic errors are removed by the procedure described above; thus, errors in the derived 
quantities, reflecting only the random noise in the input data, can be treated separately. 

(2) The errors in the derived quantities result only from errors in the Doppler residuals. The 
uncertainties in other factors entering into the analysis, such as the wavelength of the telemetry signal 
and the spacecraft position, either are negligible compared with the errors in the residuals or are 
systematic and hence removed in the data reduction. 

(3) The errors in the Doppler residuals are uncorrelated and statistically independent. 

In accordance with the notation of the previous section, if a? denotes the variance in the Doppler 
residuals at the /th level, where the superscript 5 refers to Doppler residuals, the variance in a quan- 
tity Q derived from the residuals at the nth level is 



(18) 


where Q n is the value of Q at the nth level, assumed to be a function of the Doppler residuals at level n 
and at all higher altitudes (i.e., at all levels / < n). 

For the purposes of the following discussion, it is convenient to define two matrices M and S by 


dQ; 

= , S (? 6 = 

v dS. v 



(19) 


Thus, the yth element of the matrix M^ 6 is the derivative of Q at altitude i with respect to 8 at altitude 
/, and the matrix S^ 6 is obtained by squaring the individual elements of . Since all derived quan- 
tities are obtained by downward integration from the top of the atmosphere, their values at a given 
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altitude do not depend on the Doppler residuals at lower altitudes. All of the matrices used in this 
analysis (with the single exception noted below) are, therefore, either lower triangular or diagonal. 

The specific methods used in obtaining the random error in the various atmospheric quantities derived 
from the Doppler residuals are enumerated below. 


Phase Change 


In view of Equation 2, the matrix elements involved in the calculation of the errors in the phase 
change are 


M }‘ = 


0, 

if i = 1 or i < / 

1/ ^A*/,/+ 1 ’ 

if i > 1 and / = 1 

fcCAfy-i,/ + A*/,/+ 1 ) > 

if 1 < / < i 


II 

V 


( 20 ) 


The phase-change variance is 


a* = S 06 a 6 , 


( 21 ) 


where a 0 is the vector whose components are the phase-change variances at points /, S 05 is defined by 
Equation 19, and a 5 is the (constant) vector of the residual variance. 


Straight-Line Approximation to Phase Change 


The variance in <£ is given by 


o * =S° 6 a 6 , 


( 22 ) 


where is the vector whose components are the variances of at levels /. In view of Equation 6, the 
matrix S° 6 is easily computed using Equation 19 once M^ 6 is found from the equation 


or 


3<h. 

a<t> ; 

d( t>k 

3$. 

3o * 

= 

= 

_l_ 


v 35 / 


3 §/ 

9o * 

3 5/ 


(23) 


M 06 =M O0 M 06 + M <I,a M a6 . 


(24) 


The matrices and are calculated directly from Equation 6 by using Equation 19, with 0 
and a replacing 5. We may note that = U, the unit matrix, whereas M <t>a is a constant multiple of 
the unit matrix. M 06 has already been computed in estimating the phase-change errors (Equation 20) 
and M a6 is readily obtained from Equations 2 and 4. After carrying out the matrix multiplication in 
Equation 24, we obtain 


M,r - 




Zs i a i 


D. 


(25) 


where D { . is the Kronecker delta. 


17 



Refractivity 


The variance in refractivity is given by 

o*=S Rs • o s . (26) 

The refractivities are not given explicitly as functions of the Doppler residuals but are expressed in 
Equation 1 2 as functions of the straight-line approximation to the phase change. The matrix M^ 6 , 
required for the computation of S Rs , may therefore be obtained from 

M * 6 =m ^ 4> . M a>6 9 (27) 

where is obtained from Equation 12, the recursion formulae for the refractivity at various atmos- 
pheric levels, and M 4 " 6 has been calculated in Equation 25. If we define a matrix |3 as 

X • 10 6 


and 

AZ.. 

% = "2— ^ (/>/), (28) 

the matrix * may be calculated from the equations 

M (29) 
and 



(30) 


Electron Density 


The basic equation for the variance in electron density is 

a^e = s N e s • a s . (31) 

The matrix M 7 ^ 6 , required for the calculation of S^ 6 , is obtained from 

=M N e R • M* 6 , (32) 


where is a multiple of the unit matrix computed directly from the relation in Equation 1 7 be- 
tween electron density and refractivity and M^ 5 has been calculated in Equation 27. The matrix M NeR 
is given by 


M w.e R = - 


D„ . 


40.3 X 10 


12 " 


(33) 
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Temperature 


The temperature variance is given by 

o T = S r6 • o 5 . (34) 

The matrix M 76 is obtained from 

M r 6 = M tr • M* 6 , (35) 

where M^ 6 has been calculated in Equation 27 and M TR may be computed from the expression (Equa- 
tion 16) relating temperature to refractivity. Since the temperature at the top of the lower atmosphere 
7/ 0 is computed from the scale height determined by the two uppermost refractivity points, the matrix 
M 77 * has a single nonzero element above the diagonal The matrix elements corresponding to deriva- 
tives of the temperature at the initial altitude are 


and 


mg 

M tr = Ah 

1V1 11 £ ^ n 12 



M 


TR 

12 


*1 

L Mff , 

R 2 11 


(36) 


(37) 


where g is the average value of the acceleration due to gravity in the uppermost layer A h l2 - If we use 
the trapezoidal rule to integrate the refractivities in Equation 16, the diagonal matrix elements of 
M tr are 


where 


T(h 0 )R(h 0 ) 

M TR = + 

it o 

r} 


m 





(/> 1 ), 


(38) 


A/Zj 2 

x = (g l R i +g 2 R 2 ) — r— + ■ • ■ + [g t _ x Rj-i + gjRj) — — 


(39) 


The off-diagonal elements are 




m 

2 kR t 


%- 1 Ah f- 1 ,/ + Sf A h. /+1 ) a > /) . 


(40) 


The standard deviation in a quantity at a given altitude is the square root of the variance at that 
altitude, and the 3a errors quoted below are three times the various standard deviations so obtained. 
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RESULTS 


The results presented here were obtained by processing the averaged sets of Doppler data at 
each of the four Mariner occultation points. In all four cases, analysis of the sets of data taken by the 
individual DSN stations produced results that are closely similar to those described here. 

The profiles of phase change, refractivity, electron density, temperature, and pressure shown 
were derived by considering 300 different polynomial fits to the raw residual data at each of the 
Mariner 6 and 7 occultation points. All of these fits were consistent with the zero-refractivity-ledge 
assumption, described above; i.e., they led to refractivity levels in the middle atmosphere between 
-0.05 and +0.05 «-unit. The heavy curves in Figures 1 1 through 24 correspond to the averages of 
300 profiles of the phase change and derived atmospheric parameters at each of the occultation points. 
The crosshatched areas surrounding the heavy curves in these figures correspond to the ±3a range in 
each parameter, obtained from the ensemble of the 300 cases. This range can be considered a measure 
of the systematic error produced by the imperfect removal of the nonatmospheric components in the 
raw residual data. The areas bracketing the zones of systematic error in these figures represent the 
uncertainty, in each parameter, due to the growth of formal errors in the calculations. 

The coordinates, local time, and local zenith angle at the four occultation points are given in 
Table 1. The season on Mars at the time of the Mariner 6 and 7 encounters was early fall in the 
northern hemisphere, with a solar declination of about 8°S. 

In Figures 1 1 through 14, the straight-line approximation to the phase change is shown as a func- 
tion of distance from the limb of Mars along the perpendicular to the Earth-Mars line at the four 
occultation points. In all four cases, the actual phase change differs only slightly from the straight- 
line approximation shown throughout most of the encounter period; the maximum difference is about 
30 percent at the surface, in the case of Mariner 6 exit. The systematic and formal errors in the actual 
phase change are also closely similar to those bracketing the curves shown. This is a result of the small 
bending angle produced by the thin Martian atmosphere. This angle reaches a maximum of about 
0.35 mrad at the surface of Mars, in the case of Mariner 6 exit. No evidence of an ionospheric-induced 
phase change was found in the Mariner 6 and 7 exit cases, which occurred at nighttime on Mars; and 
in these cases, only the phase change due to the lower atmosphere is shown (Figures 1 2 and 1 4). The 
absence of a region of detectable negative refractivity at the exit points in the upper atmosphere of 
Mars results in a larger positive phase change as occultation is approached. The presence of a spherical 
shell of negative refractivity at high altitudes with the properties predicted for the Martian ionosphere 
(Stewart and Hogan, 1969) would contribute negatively by about 10 cycles to the phase change 
experienced by a ray passing near the surface in a one-way transit of the atmosphere. 

Profiles of atmospheric refractivity at the four occultation points are shown in Figures 1 5 
through 18. In each case, these curves are averages of the 300 refractivity profiles obtained by inver- 
sion of the 300 individual phase curves whose average properties are depicted in Figures 1 1 through 
14. The error zones surrounding the refractivity profiles were obtained as discussed above. 

Averages of 300 profiles of electron density in each of the Mariner 6 and 7 entry cases are pre- 
sented in Figures 1 9 and 20, along with their accompanying systematic and formal (± 3a) error limits. 
As can be seen, the uncertainty in electron density resulting from formal errors is larger by an order 


20 



PHASE CHANGE (cycles) x ^ — PHASE CHANGE (cycles) 



gure 11— Phase change vs. distance from limb of Mars for 
iriner 6 entry case. Error bars are ±3 a. (See text for 
planation.) 
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Figure 12— Phase change vs. distance from limb of Mars for 
Mariner 6 exit case. Error bars are ±3a. (See text for 
explanation.) 




Figure 13— Phase change vs. distance from limb of Mars for 
Mariner 7 entry case. Error bars are ±3a. (See text for 
explanation.) 



Figure 14— Phase change vs. distance from limb of Mars for 
Mariner 7 exit case. Error bars are ±3a. (See text for 
explanation.) 









ELECTRON NUMBER DENSITY 

Figure 19— Electron number density (cm -3 ) vs. altitude for 
Mariner 6 entry. Error bars are ±3a. (See text for 
explanation.) 
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Figure 20— Electron number density (cm -3 ) vs. altitude for 
Mariner 7 entry. Error bars are ±3a. (See text for 
explanation.) 



Table 1 — Occultation parameters. 






Local Solar 

Occultation 

Latitude 

Longitude 

Local Time 

Zenith Angle 





(deg) 

Mariner 6 entry 

3.7°N 

355. 7°E 

1545 

57 

Mariner 6 exit 

79.3°N 

87.1°E 

2210 

107 

Mariner 7 entry 

58.2°N 

30.3°E 

1430 

56 

Mariner 7 exit 

38.1°N 

21 1.7°E 

0310 

130 


of magnitude than that resulting from systematic errors. For the Mariner 6 and 7 entry cases, the 
average peak electron densities are 1.59 cm -3 at 136.5 km and 1.73 cm -3 at 137.5 km, respectively. 
The formal error limits (±3a) on the electron density at the peak are about ±2.0 X 10 4 cm -3 in both 
cases. Determinations of the scale height of ionization between 150 and 225 km in these cases were 
made by fitting a straight line to the logarithms of the electron densities at intermediate levels. For 
the average of 300 cases, the topside scale heights are 43.0 km and 47.0 km, respectively, with formal 
(±3a) uncertainties of ±6.0 km and ±4.0 km. With the assumption of a 100 percent C0 2 atmosphere 
at these levels, the corresponding topside temperatures are 388 ± 54 K in the case of Mariner 6 entry, 
and 424 ± 35 K in the case of Mariner 7 entry, where the error bars on temperature correspond to the 
formal (±3a) errors. The larger error limits on the topside temperature in the Mariner 6 entry case are 
due to the presence of a dip in the electron profile above about 200 km (Figure 1 9). 

In Figures 21 through 24, the averages of 300 temperature profiles and the attendant errors are 
shown for each of the four occultation points, for an assumed atmospheric composition of 1 00 per- 
cent C0 2 • In these figures, the profile of the saturation temperature of C0 2 is shown as circles. In 
Table 2, the values of atmospheric temperature and pressure at the surface of Mars, with the estimated 
uncertainties due to both systematic and formal errors, are given for the four occultation points. 
Ionospheric parameters are also indicated in this table. 

The effect of the adopted values of the occultation times on the derived values of the atmospheric 
parameters at the surface was estimated by varying this quantity by ±0. 1 s. In all cases, the tempera- 
ture at the surface varied by less than ±2 K, and the pressure at the surface varied by less than ±0.5 mb. 

Temperature distributions were also derived for 95 percent C0 2 , 5 percent N 2 and 90 percent 
C0 2 , 10 percent Ne mixtures. The results of these calculations, also averaged over 300 cases, are 
shown in Figures 25 through 28, along with the average profiles obtained for the 100 percent C0 2 
cases. In deriving these profiles, the temperature was constrained by the saturation temperature of 
C0 2 . Since the observed refractivity scale height essentially fixes the ratio of temperature to mean 
molecular mass, the effect of the presence of lighter gases such as N 2 or Ne in the atmosphere is to 
lower the temperature while increasing the pressure. For clarity, the error zones have been omitted 
from these figures. The errors are of magnitude similar to the corresponding cases shown in Figures 
21 through 24. 
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MARINER 6 ENTRY 
AVERAGE CASE 



Figure 21— Temperature vs. altitude at Mariner 6 entry. 
Saturation temperature of C0 2 is shown as circles. 
Systematic error bars are ±3a; formal error bars are ±o. (See 
text for explanation.) 



Figure 22— Temperature vs. altitude at Mariner 6 exit. 
Saturation temperature of C0 2 is shown as circles. 
Systematic error bars are ±3a; formal error bars are ±o. (See 
text for explanation.) 
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Figure 23— Temperature vs. altitude at Mariner 7 entry. 
Saturation temperature of C0 2 is shown as circles. 
Systematic error bars are ±3a; formal error bars are ±o. (See 
text for explanation.) 
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Figure 24-Temperature vs. altitude at Manner 7 exit. 
Saturation temperature of C0 2 is shown as circles. 
Systematic error bars are ±3a; formal error bars are ±o. (See 
text for explanation.) 




Table 2— Atmospheric parameters on Mars from Mariners 6 and 7. 


Event 

Surface 

Temperature* 

(K) 

Surface 

Pressure* 

(mb) 

Maximum 
Electrons** 
(X 10 5 cm -3 ) 

Peak 

Altitude 

(km) 

Topside 

Scale 

Height** 

(km) 

Topside 

Temperature** 

(K) 

Mariner 6 entry 

254 ±7 

6.07 ± 0.36 

1.59 ±0.2 

136.5 

43± 6 

388 ± 54 

Mariner 6 exit 

152 ± 5 

8.51 ±0.48 

— 

— 

— 


Mariner 7 entry 

221 ± 8 

4.95 ± 0.24 

1.73 ±0.2 

137.5 

47 ±4 

424 ± 35 

Mariner 7jexit 

— 

209 ±6 

8.02 ± 0.45 

— 

— 


— 


*100 percent C0 2 atmosphere. Error bars correspond to ±o formal errors plus ±3 o systematic errors. 
**Error bars correspond to ±3o formal errors plus ±3 a systematic errors. 



ALTITUDE (km) 

Figure 25— Temperature vs. altitude for Mariner 6 entry. 
Saturation temperature of C0 2 is shown as circles. Curve 
1—100 percent C0 2 ; Curve 2—95 percent C0 2 , 5 
percent N 2 ; Curve 3—90 percent C0 2/ 10 percent Ne. 
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Figure 26— Temperature vs. altitude for Mariner 6 exit. 
Saturation temperature of C0 2 is shown as circles. Curve 
1—100 percent C0 2 ; Curve 2—95 percent C0 2 , 5 
percent N 2 ; Curve 3—90 percent C0 2 , 10 percent Ne. 
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Figure 27— Temperature vs. altitude for Mariner 7 
entry. Saturation temperature of C0 2 is shown as 
circles. Curve 1—100 percent C0 2 ; Curve 2—95 per- 
cent C0 2 , 5 percent N 2 ; Curve 3—90 percent C0 2 , 
10 percent Ne. 


Figure 28— Temperature vs. altitude for Mariner 7 
exit. Saturation temperature of C0 2 is shown as 
circles. Curve 1 — 100 percent C0 2 ; Curve 2—95 per- 
cent C0 2 , 5 percent N 2 ; Curve 3—90 percent C0 2 , 
10 percent Ne. 
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The pressure distributions at the four occultation points (averages of 300 cases) for an atmosphere 
of 100 percent C0 2 are shown in Figure 29. 

DISCUSSION 

The Mariner 6 and 7 temperature profiles shown in Figures 21 through 24 are similar to those 
published by Rasool et al. (1970). The principal difference is in the Mariner 6 exit profile, which is 
now about 10 K colder at all altitudes than the corresponding profile obtained in the earlier analysis. 

As noted above, the results of Rasool et al. were based on a single reduction of the residual data and 
employed a different bias-removal technique. Since the present results are based on more extensive 
analyses of the data, they should be regarded as more reliable. 

The temperature profiles of Figures 21 through 24 contain several features of interest. In all of 
the thermal profiles derived from the Mariner 6 and 7 occultation experiments, temperatures in the 
middle atmosphere fall below the value at which C0 2 would saturate. This drop in temperature is 
sharp; in all four cases, lapse rates in excess of adiabatic are found throughout an atmospheric layer 
several kilometers thick. In some cases, lapse rates as high as 10 K/km are obtained. 

The indication of C0 2 saturation in the middle atmosphere is the most pervasive single feature of 
the derived temperature profiles; it occurs regardless of the details of the raw data adjustment. If the 
occultation data alone suggested C0 2 saturation in the middle atmosphere, this result would have to be 
viewed with some skepticism; the relative errors in refractivity are largest in this region where refractiv- 
ities are small and the corresponding errors in the derived temperatures are large (Figures 21 through 
24). The fact that all of the temperature profiles show this effect would imply that it is real, although 
an alternate explanation, in terms of the possible existence of a Martian D-region ionosphere, was 
pointed out by Rasool et al. (1970). There are other indications, however, apart from the occultation 
data, which suggest that C0 2 condensation does occur in the middle atmosphere. 

The infrared spectrometers aboard both Mariners 6 and 7 recorded a reflection at 4.3 jum from 
altitudes well above the planetary surface (Herr and Pimentel, 1970). At the second Mariner 7 bright- 
limb crossing, the experimenters estimate the altitude of the reflecting feature to be 25 ± 7 km. This 
reflection spike indicates the presence of condensed C0 2 in the Mars middle atmosphere. 

A third line of evidence suggesting possible C0 2 condensation comes from the Mariner 6 and 7 
TV pictures (Leovy et al., 1971), where thin haze layers at altitudes ranging from 5 to 50 km above 
the limb of Mars were observed in several frames. Some of these layers are high enough in the atmos- 
phere so that the scattering could be due to solid C0 2 , although, as noted by Leovy et al. (1971 ), the 
lower layers seem to require another mechanism. 

As pointed out above, prior to the Mariner 6 and 7 experiments, no theoretical models of the 
Martian lower atmosphere predicted condensation of C0 2 . However, Gierasch (1971 ) has recently 
suggested that the kinetic energy generated by buoyancy might not be dissipated locally on Mars; this 
would result in the penetration of turbulent motions into what would otherwise be a stable radiative 
regime. Such overshooting would raise the tropopause level, producing substantially lower tempera- 
tures in the middle atmosphere. According to Gierasch (1971), at low latitudes the temperature would 
be driven low enough by this process for C0 2 to condense. The present results indicate, however, that 


31 



PRESSURE (mb) 


10-2 b" 


1— Mariner 6 Entry 

2— Mariner 6 Exit 

3— Mariner 7 Entry 

4— Mariner 7 Exit 


0 10 20 30 40 5 

ALTITUDE (km) 

Figure 29-Pressure vs. altitude at the four occultation points for an atmosphere of 100 percent C0 2 . Range in observed surface 
pressure levels corresponds to about a 4-km altitude range, based on a pressure scale height of 10 km. 



condensation occurs at high as well as low latitudes, and that the temperature gradient in the lowest 
20 km of the atmosphere is subadiabatic. Neither of these findings is consistent with the theory of 
Gierasch, and it appears that at present there is no complete explanation for the thermal structures 
indicated in Figures 21 through 24. 

The temperature obtained here for the atmosphere near the surface at the Mariner 6 entry point 
is 254 ±7K, whereas the ground temperature obtained with the infrared radiometer for the same 
region is 275 ± 5 K (Neugebauer et al., 1969 and 1971). Thus, an air-ground temperature discontinuity 
of about 20 K in the late Martian afternoon is indicated. The existence of a sizable temperature dis- 
continuity between the atmosphere and surface of Mars in the afternoon was predicted by Gierasch 
and Goody (1968). According to these authors, turbulent transfer is the dominant process for heat 
exchange in the lower atmosphere. For the estimates of the turbulent heat flux which they used, 
Gierasch and Goody obtained a maximum temperature discontinuity (ground minus air) of +65 K for 
early afternoon and -40 K for early morning hours. The present results are in agreement with the 
prediction that a sizable temperature discontinuity can exist at times at the surface of Mars. However, 
the observed discontinuity of +20 K, in contrast with the discontinuity of +40 K predicted for 1600 
local time, indicates that the magnitude of turbulent heat transfer in the Martian boundary layer is 
probably higher than Gierasch and Goody assumed. 

In the same context, the Mariner 6 exit profile, corresponding to 2200 local time, exhibits a 
strong inversion near the ground (Figure 22), with the maximum temperature occurring near 4 km, 
about 15 K higher than the temperature at the surface. A similar feature was obtained by Rasool 
et al. (1970) but was not present in the profiles of Kliore et al. (1970). On the other hand, the night- 
side profile at Mariner 7 exit does not exhibit an inversion (Figure 24), and, according to the time- 
dependent calculations of Gierasch and Goody (1968), all of the nightside profiles should be charac- 
terized by such a feature. However, the absence of an inversion in the Mariner 7 exit profile cannot 
be interpreted with certainty as a disagreement with theory, since the error bars derived for this case 
permit a mild inversion. 

The observed differences in total pressure at various points on Mars (Table 2) are almost certainly 
due to the topography of the surface; the lower pressure values correspond to elevated regions, where- 
as the higher pressures indicate depressed areas. Because the pressure scale height in the lower atmos- 
phere of Mars is about 10 km, the observed pressure extremes would correspond to altitude variations 
on the surface of as much as 4 km. These variations are consistent with the recent radar results of 
Pettengill et al. (1969) and with the variation in C0 2 abundance over the disk of Mars observed both 
in ground-based measurements (Belton and Hunten, 1971) and in the Mariner 6 and 7 data (Herr et al., 
1970). 

The basic results of the Mariner 6 and 7 occultation experiments for the upper atmosphere of 
Mars are the profiles of the distributions of electron density with altitude (Figures 19 and 20). Several 
points about these ionospheric profiles should be noted. The maximum electron density observed at 
the time of the Mariner 4 flight (near solar minimum) was (9.0 ± 1.0) X 10 4 cm” 3 ; this was found 
near 1 20 km (Fjeldbo and Eshleman, 1968). Four years later, during a period of intermediate solar 
activity, Mariners 6 and 7 observed electron densities between (1.6 ± 0.2) X 10 5 cm* 3 and 
(1.7 ± 0.2) X 10 5 cm""" 3 near 135 km, an increase in ionization by nearly a factor of 2. Similar values 
were derived by Fjeldbo et al. (1970) for Mariners 6 and 7. 
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A factor to be noted here is the departure of the Martian ionosphere from spherical symmetry. 
Such a departure would be most important in the interpretation of the Mariner 6 data because the 
occulting ray in this case passed across the occultation point in the general direction of the subsolar 
point. As the subsolar point is approached, the locus of the maximum electron density decreases in 
altitude, whereas the magnitude of the local maximum increases. The total number of electrons along 
the ray path would thus be greater for rays passing below the maximum electron density at the occul- 
tation point. Thus, the altitude of the peak ionization obtained by assuming spherical symmetry in 
the inversion of the refractivity data would be below the actual altitude of peak ionization at the 
occultation point. On the basis of numerical experiments in which Chapman-type ionospheres with 
different peak electron densities, scale heights, and peak altitudes were fitted to the phase change data, 
the authors estimate that the error in the observed altitude of maximum ionization due to this effect 
would be between 1 and 3 km. 

The primary electron-density maxima observed in the Mars upper atmosphere are now generally 
accepted to be Fj ionospheric layers, although there is still some controversy as to whether they are 
in strict photochemical equilibrium (R. W. Stewart, 1971) or are modified via interaction with the solar 
wind (Cloutier et ah, 1969). If the ultraviolet fluxes reported by Hall and Hinteregger (1970) are 
assumed to be the only source of ionization, the observed electron-density profiles exhibit excess 
ionization by as much as a factor of 2 over theoretical models (R. W. Stewart, 1971 ; Donahue, 1971). 

If the major electron-density maximum in the Mars ionosphere is interpreted as an F 1 -layer in 
photochemical equilibrium, the plasma scale heights above the peak can be used to infer an exospheric 
temperature. At the time of the Mariner 4 flight, the exospheric temperature obtained in this way was 
about 275 K (Fjeldbo and Eshleman, 1968). The Mariner 6 and 7 plasma scale heights of 43 km and 
47 km, measured between altitudes of about 150 and 225 km at a time of increased solar activity, 
yield exospheric temperatures of 388 ± 54 K and 424 ± 35 K, respectively. Similar results were found 
by Fjeldbo et al. (1970). Also, there is an indication of a large increase in the plasma scale height 
above 225 km at both the Mariner 6 and 7 entry points. Such an increase was also found by Fjeldbo 
et al. (1970), who suggested that it might represent a change of the dominant ion species in an iono- 
sphere in diffusive equilibrium. It may also be attributed to a transition from an ionospheric region in 
photochemical equilibrium to a region where ambipolar diffusion controls the distribution of ions. 

Both these plasma scale heights (below 225 km) and the neutral scale heights of about 25 km 
deduced from the Mariner ultraviolet spectrometer data (A. I. Stewart, 1972) indicate substantially 
lower exospheric temperatures than were predicted on theoretical grounds by McElroy (1969) and 
Stewart and Hogan (1969). Also, the fact that the derived plasma scale height is approximately twice 
the neutral scale height observed implies the absence of any alteration of the plasma distribution by 
the solar wind. 

In addition to the major ionization peak, all of the Mariner flights to Mars (and Venus) have 
observed a secondary electron ledge at lower altitudes (Figures 19 and 20). This ledge appears about 
25 km below the primary maximum in the Mariner 6 and 7 cases and was observed in a similar position 
by Mariner 4. This secondary electron peak is analogous to the Earth’s E-region. It appears naturally 
in theoretical models of the Mars and Venus ionospheres (e.g., R. W. Stewart, 1971) and is a con- 
sequence of ionization of atmospheric constituents by soft X-rays. 
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